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ABSTRACT 

In  a  reacting  gas  mixture  the  frozen  sound  speed,  a  ,  is 

epfp.s.y.]  1 

a2  =  (-«^ - )S,  y.  while  the  equilibrium  sound  speed,  a^,  is  given  by 

dP[p  ,  S,  yie(p.  Sfl 

a|  =  ( - ^ - -)g.  The  expressions  for  the  two  sound  speeds  are  not  equal, 

indicating  a  discontinuity  as  the  reaction  rate  approaches  infinity.  In  this  report 
the  one-dimensional  propagation  of  sound  and  the  significance  of  the  two  sound 
speeds  is  examined  in  detail  for  the  case  of  a  gas  mixture  at  or  near  equilibrium, 
consisting  of  gases  which  obey  the  ideal  gas  law  and  Dalton' s  law  of  partial  pres¬ 
sures  and  which  have  all  their  degrees  of  freedom  fully  excited.  The  significant 
results  obtained  are  summarized  below: 


1)  The  movement  of  a  pressure  perturbation,  5P,  is  characterized  by 


1  8^(6P)  92(6P) 
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to  explain  the  physical  meaning  of  the  two  sound  speeds,  a  and  af. 

2)  For  an  instantaneous  reaction,  the  state  and  composition  variables  are 
shown  to  be  representable  as  solutions  of  the  characteristic  differential  equations 
and  as  step  functions,  and  the  rate  of  reaction  as  an  impulse  function. 

3)  Explicit  expressions  are  obtained  for  ag  and  a^  in  terms  of  the 
equilibrium  state  and  composition  variables. 

NOMENCLATURE 

a  speed  of  sound  of  nonreacting  medium  = 

dPp  ,  S,  y  (p,  S)  2 

a^  equilibrium  speed  of  sound  =  (  ■  - . —  )  g 

a^  frozen  speed  of  sound  = 

any  possible  component  of  the  gas  mixture 
C.  undetermined  coefficient 

l 

C  constant  pressure  specific  heat  of  gas  mixture,  energy  per  unit  mass  per  < 

P 

temperature  degree 

C  constant  pressure  specific  heat  of  A 

C  constant  volume  specific  heat  of  gas  mixture 

V  op 

e  base  of  natural  logarithms  =  )  -r-j-  «  2.  72 

i-1  1 ' 

HL  specific  enthalpy  of  A^,  which  includes  heat  of  formation 
AHr  enthalpy  change  for  reaction 

J  mechanical  equivalent  of  heat 

k  specific  heat  ratio 

Kc  equilibrium  constant  in  terms  of  molal  concentrations 
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equilibrium  constant  in  terms  of  partial  pressures 
average  molecular  weight  of  gas  mixture 
molecular  weight  of  component 

molal  density  of  the  gas  mixture,  total  number  of  mols  per  unit  volume 

pressure  of  gas  mixture  at  any  point 

speed  of  gas  stream 

universal  gas  constant 

specific  entropy  of  gas  mixture 

specific  entropy  of  at  T  and  P 

specific  entropy  of  A^  at  any  T  and  P 

entropy  change  for  reaction 

time 

iefined  as  t  =  t 
r 

absolute  temperature 

mass  rate  of  production  of  A.  =  y 

distance  coordinate  along  the  direction  of  flow 

defined  as  x  =  x  -  at 
r  ^ 

mass  fraction  of  A^  in  the  gas  mixture 
, S)  mass  fraction  of  A^  under  equilibrium  conditions 
total  number  of  possible  constituents  in  the  gas. mixture 
undetermined  coefficients 
stoichiometric  coefficient  of  A..  in  reaction 
net  rate  of  reaction  defined  in  Eq.  (12) 
parameter 
determinant 

free  energy  of  component  A^,  defined  as  =  FL  -  TS^ 
free  energy  change  of  reaction 
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p  mass  density  of  gas  mixture 

T  time  of  reaction 

backward  rate  constant  in  Eq.  (32) 

T £  forward  rate  constant  in  Eq.  (32) 

y- 

T*  defined  by  Eq.  (33) 

T  defined  by  Eq.  (65) 

>  means  greater  than 

<  means  less  than 

Functional  and  Operational  Symbols 
differential  of  (  ) 
perturbation  of  (  ) 

partial  derivative  of  (  )  with  respect  to  [  ] 

unit  step  function  defined  by  U(\)  =  0  for  \<  0  and  U(\)  =  1  for  \>0 

00 

impulse  function  defined  by  I(\)  =  0  for  \  /=  0  and  \  I(\)d\  =  I 

J_oo 

function  defined  by  <\>=  \  for  \>  0  and  <\>  =  0  for  \  <  0 
integration  sign 
summation  sign 
product  sign 

usually  runs  from  1  to  Z 
Superscript 

—  denotes  quantity  at  equilibrium 
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INTRODUCTION 


Dr.  Boa-Teh  Chu  has  pointed  out  that  the  expression  for  the  speed  of 
sound  in  a  reacting  gas  mixture  appears  to  have  a  discontinuity  as  the  reaction 
rate  approaches  infinity.  By  setting  up  differential  equations  for  the  gas  flow,  it 
is  possible  to  show  by  the  method  of  characteristics  that  the  sound  speed,  a^,  (com¬ 
monly  known  as  the  frozen  speed  of  sound)  is  given  by 


af3  ■  < 


3P  [p,S,y.] 


where  P  is  the  total  pressure,  pis  the  density,  S  is  the  specific  entropy  of  the 
gas  mixture,  and  y.  (i  =  1,  2,  3, . . . . ,  Z)  is  a  composition  variable.  As  indicated, 
the  differentiation  is  carried  out  holding  the  entropy  and  composition  variables 
constant.  This  expression  is  in  agreement  with  that  for  a  nonreacting  medium, 

af  = 

It  is  known,  however,  that  for  a  gas  mixture  undergoing  a  very  fast  reac¬ 
tion  (i.  e.  ,  one  in  which  the  time  of  reaction  approaches  zero)  the  sound  speed,  a  , 
(known  as  the  equilibrium  speed  of  sound)  is  given  by 

3P[p,S,y  •  (p,S)] 

e  '  (  >S  (3) 

where  y.  (o,S)  is  the  equilibrium  composition,  which  depends  on  o  andS,  and  the 
differentiation  is  carried  out  holding  only  the  entropy  constant.  The  two  expressions 
are,  in  general,  not  equal  and  are,  in  fact,  related  by 


a2  _a3  +  V 

ae  '  a£  +.L, 

1=  1  'l 


3p  8yie<°’s> 
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where  Z  is  the  number  of  constituents  in  the  mixture.  Moreover,  since  the  reac¬ 
tion  rate  does  not  appear  in  the  expressions  for  a^  and  a^,  there  is  no  way  to  obtain 
one  from  the  other  by  a  limiting  process. 

In  this  paper  the  phenomenon  of  the  two  sound  speeds  is  examined  in  detail 
for  the  case  of  one -dimensional  gas  flow  in  which  the  gas  mixture  is  near  chemical 
equilibrium  and  in  which  it  is  assumed  that  all  the  degrees  of  freedom  of  each 
molecule  are  fully  excited.  This  last  condition  is  closely  approximated  at  high 
temperatures  and  high  speed  gas  flow. 

A  general  set  of  differential  equations  is  obtained  for  the  gas  mixture  near 
equilibrium  and  then  the  effect  of  a  reaction  with  a  zero  and  a  non- zero  reaction 
time  is  studied.  It  is  shown  that  the  characteristic  surfaces  of  the  differential 
equations  move  at  speed  a^  for  all  non-zero  reaction  times,  but  that  the  speed  of 
the  characteristics  changes  dis continuously  to  a  in  the  (physically  unrealizable) 
case  of  a  reaction  time  equaling  zero. 

A  differential  equation  is  then  derived  for  the  propagation  of  a  pressure 
pulse  through  the  gas  stream,  which  clarifies  the  relationships  between  a  ,  a,, 
and  the  time  of  reaction.  In  particular,  it  is  deduced  from  the  above  equation  that 
for  a  small  reaction  time,  a^  represents  the  speed  of  the  wave  front,  while  ag 
represents  the  velocity  of  the  bulk  of  the  gas  behind  the  wave  front,  which  is  in 
agreement  with  the  results  given  in  reference  1. 

The  theoretical  problem  of  an  instantaneous  reaction  is  discussed  and  it  is 
shown  that  this  problem  can  be  solved  by  considering  the  state  and  composition 
variables  to  be  step  functions  and  the  reaction  rate  to  be  impulse  function. 

In  the  course  of  this  work  explicit  expressions  are  obtained  for  a  and  a  . 


GENERAL  THERMODYNAMIC  REALTIONSHIPS 


Suppose  Z  different  gases  are  moving  in  one -dimensional  flow  through  an 
insulated  duct  of  constant  cross-section  with  a  velocity  q.  In  general,  several 
reactions  can  occur  between  some  or  all  of  these  gases,  the  predominant  reactions 
depending  on  the  temperature,  pressure,  and  composition  at  any  instant.  Even¬ 
tually,  the  gases  will  reach  an  equilibrium  state  and  the  composition  and  internal 
state  variables  for  this  equilibrium  state  can  be  calculated  from  thermodynamic 
considerations  and  the  initial  state  of  the  gases.  If  this  equilibrium  is  disturbed 
slightly,  for  example,  by  a  small  pressure  or  temperature  change,  an  overall 
reaction  occurs,  which  can  be  written  as 


z 


y.a.  =  0 
1  1 


(5) 


where  A.  is  a  constituent  and  the  stoichiometric  coefficient  of  this  constituent, 
l  ‘i 

The  y^‘s  can  ta^e  on  anY  real  value,  including  zero  and  negative  values.  If  y^  is 

positive,  the  component  Ai  is  called  a  product;  if  y^  is  negative,  A^  is  a  reactant; 

and  if  •  .  is  zero.  A.  is  a.  nonreactant  or  dilutant.  Since  the  reaction  is  reversible, 
■ii 

the  distinction  between  a  reactant  and  a  product  is  arbitrary. 

At  any  temperature  and  pressure,  component  A^  possesses  internal  thermo¬ 
dynamic  state  properties,  H.,  p. ,  where  R  is  specific  enthalpy,  S .  is  specific 
entropy,  and  p^  is  specific  free  energy.  These  properties  are  related  by 


>*i  ■  Hi  -  TSi 


(6) 


where  T  is  the  absolute  temperature.  In  general,  R,  S^,  and  p^  are  functions  of 
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temperature  and  pressure,  but  for  gases  obeying  the  ideal  gas  law,  HL  is  a  function 
of  temperature  only. 

At  any  fixed  absolute  temperature  T  and  pressure  P,  the  stoichiometric 
reaction,  Eq.  (5),  possesses  enthalpy,  entropy,  and  free  energy  changes,  AH^AS^, 
and  AUr>>  which  are  related  by 

XV 

A|iR  =  AHr  -  TASR.  (7) 


The  enthalpy  change  of  the  reaction  is  related  to  the  specific  enthalpies  of  the 
constituents  by 

z 

4Hr  =  Y.m.H.  '  (8) 

i=l 

where  rm  is  the  molecular  weight  of  A^.  Similarly, 


and 


z 

\ 


L 

i=l 


.m.Li. 
i  li 


aSr 


\ 

L 

i=l 


Y  .m.S.  . 

ill 


It  is  a  property  of  the  free  energy  that  under  equilibrium  conditions, 


(9) 


(10) 


APr 


\ 

L 

i=l 


Y  -m.Ll. 
T  l  ri 


=  0 


(Ha) 


A  Hr  =  T  A 5r  . 


(lib) 


and  so 


(From  this  point  on,  equilibrium  quantities  are  represented  by  a  bar. ) 
The  rate  of  reaction,  p  >  is  defined  as 


r  = 


i  d'Ai> 


v.  dt 

’i 


(12) 


d(A. ) 


1 

where  — —  represents  the  molal  rate  of  production  of  A^.  The  mass  rate  of 
production  of  a  gas  A^  is  called  and  is  given  by 


w.  -  v.m.r  . 


(13) 


Necessarily, 


Z i 

1-,- 


(14a) 


i=l 


and  so 


A 

/  y  .m.  =  0  . 
/  l  i 


(14b) 


i=i 


Eqs.  (14)  merely  assert  that  the  stoichiometric  Eq.  (5)  is  balanced. 


GENERAL  DIFFERENTIAL  EQUATIONS  OF  THE  SYSTEM 

The  mass  balance  in  a  gas  stream  for  any  one  of  the  components  is  written  as 


8(py;)  8<pyi>  8q 

_ _  n _ 4-  \r _ — 


9t 


+  q 


9x 


+  P^i 


=  w. 


9x 


(15) 


where  p  is  che  mass  density  of  the  gas,  the  mass  fraction  of  A^  in  the  gas 
mixture,  and  q  the  speed  of  the  gas  stream.  Summing  Eq,  (15)  from  1  to  Z  and 


z 

remembering  that  y^  =  1,  and 

i=l  i=l 


z 

^  w.  =  0,  gi 


gives 
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3p  ,  3p  8q  A 

"5t  +  q  S  +  "  °* 


If  the  product  differentiation  in  Eq.  (15)  is  carried  out  explicitly  and  combined  with 


Eqs.  (13)  and  (16),  then 


wi  Yimir 


w + 


A  momentum,  balance  gives 


1  3P 
p 


3q  ,  Q  3q 
"SF  +  q  ‘5x 


where  P  is  the  pressure  of  the  gas  mixture  at  any  point. 

The  equation  of  state  for  the  gas  mixture,  assuming  an  ideal  gas  and 


Dalton's  Law,  is 


z 

\  yi 

P  =  nR  T  =  R  Tp  >  — 
o  o  H  /  m. 

i=l  1 


where  Rq  is  the  universal  gas  constant,  T  is  the  absolute  temperature  of  the 
system,  and  n  is  the  molal  density,  i.  e.  ,  the  total  number  of  mols  of  gas  per  unit 
volume. 

Choosing  the  equilibrium  temperature,  T,  and  pressure,  P,  as  reference 
states  at  which  each  component  has  an  entropy  per  unit  mass  of  S^q,  the  entropy 
per  unit  mass  of  gas  mixture  at  any  temperature  and  pressure  close  to  T  and  P 
is  given  by 


\  f  „  R  _  m. 

S  =  /  y .  (s .  +  C  .  In  —  -  — ^  In  — — - — 

L  M  »  pl  T  mi  p  VZ 


H 
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assuming  again  the  ideal  gas  and  Dalton's  Law. 


m. 


The  term  -  represents  the  partial  pressure  of  A.  in  the  mixture. 

.z  Y,-  1 


l 

i=l 


m. 

1 


C  .  ,  the  constant  pressure  specific  heat  of  at  equilibrium,  is,  in  general, only 
a  weak  function  of  pressure  and  temperature  and  for  a  small  pressure-temperature 
range  can  be  taken  as  constant. 

An  energy  balance  yields 


as  ,  as 
aF  +  = 


(21) 


Eqs.  (16)-(21)  are  regarded  as  six  partial  differential  equations  with  seven 
unknown  dependent  variables  y^  SiQ,P,  q,  T,  S,  and  r  and  two  independent  vari¬ 
ables  x  and  t.  is  considered  to  be  a  known  function  of  temperature  and  pressure 
and  w.  is  given  by  Eq.  (13).  In  order  to  obtain  independent  solutions  for  each  of  the 
seven  dependent  variables,  another  equation  is  needed.  This  other  equation  is  the 
relationship  between  T  >  p,  T,  and  y^  which  is  given  below. 


LINEARIZATION  AND  SIMPLIFICATION  OF  THE  DIFFERENTIAL  EQUATION 

The  differential  equations  of  the  system  are  highly  nonlinear  but  if  attention 
is  confined  to  a  section  of  the  stream  in  which  the  gas  mixture  is  very  close  to  or  at 
equilibrium,  and  if  all  disturbances  within  this  section  are  of  small  amplitude,  the 
equations  can  be  linearized  by  the  method  of  perturbations.  Using  this  method, 
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P(x,  t)  =  V  +  §P(x,  t)  y^x,  t)  =  y.  +  6Yi(x,  t) 

p(x,  t)  =  ~p  +  §p(x>  t)  T(x,  t)  =  T  +  &T(x,  t) 

q(x,  t)  =  q”  +  6q(x,  t)  S(x,  t)  =  "S'  4  6S(x,  t) 

and  (i^T.P)  =  fl.  +  6^(T,P)  . 

The  barred  quantities  are  constants  and  represent  the  values  at  thermodynamic 
equilibrium.  The  )  terms  are  perturbations  and  are  considered  to  be  of  a  much 
smaller  order  than  the  equilibrium  quantities.  Note  that  by  the  definitions  of  w.. 
and  r>  w.  and  p  are  both  zero. 

The  expressions  above  are  now  introduced  into  Eqs.  (16)-(21).  Perturbation 
terms  are  neglected  in  comparison  with  equilibrium  terms,  and  the  following  rela¬ 
tionships  are  used  for  simplification: 

w.  =  v.m.r  (22) 

1  ‘xi 

z  y. 

_  -  \  L± 

n  =  o  ^  m.  (23) 

i=l 

z 

C  =  )  C  .y.  (24) 

P  /_  pi'i 

i=l 


(25) 

(26) 

(27) 

(28) 


where  n  is  the  molal  density  at  equilibrium,  £T  the  constant  pressure  specific  heat 

P 

for  the  gas  mixture  at  equilibrium,  m  the  average  molecular  weight  of  the  gas 
mixture  at  equilibrium,  and  7>  the  specific  entropy  of  A.  in  the  gas  mixture  at 
equilibrium. 

The  linearised  equations  become 

^  +  *  ^  +  P  tgfi-  =  0  U9a) 


3(6Yi)  _  8(6^)  Y^. 

-  +  q  -  =  — — . r 

8t  0x  n 


(29b) 


3(6q)  x  —  3(6q)  A  1  8(6P)  _  A 

5t~  +  - ~ 0 


(29c) 


ifi. 

P  P 


vz  6yi 

+  il+m)  HT 

T  A.  1 

i=l 


(294) 


8(fiS)  x  -  8(6S)  _ 
“ W  +  q  53T  " 


(29e) 


(29f) 


T  is  of  the  order  of  the  derivative  of  a  perturbation  by  (29b).  Then,  since 
the  expression  on  the  left  of  (29f)  contains  derivatives  of  6,  while  the  expression 
on  the  right  contains  the  product  of  f  and  perturbation  terms,  the  expression  on 
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the  right  is  of  a  higher  order  and  can  be  neglected,  giving 


3(6S)  -8(fiS)  _ 

~W~  +  "  0 


(29g) 


The  above  equations  can  now  be  simplified  somewhat  by  changing  to  x^,  t 


coordinates  where 


x  =  x  -  qt 
r 


(30a) 


t  =  t 
r 


(30b) 


Physically,  x^  represents  the  coordinate  of  a  pulse  relative  to  a  fixed  particle  in 


the  fluid. 


Eqs.  (29)  now  become 


3(6p)  ,  -  3(6q)  _  n 

“5t“  +  p  - U 

r  r 


3(6yi)  Yjm. 


(31a) 


(31b) 


3( 6q)  A  1  0(6P)  _  n 

t  “ttsf —  u 

r  y 


(31c) 


+il  +  m  } 

p  p  t  4  mi 


(31d) 


6S  =  )  S.  6y,  +  ~  6T  -  ^  6P 
A  T  p  T 


(31e) 


8(6S) 


(31f) 
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The  derivative  with  respect  to  t  is  actually  a  "substantial"  time  derivative  and 
represents  differentiation  with  respect  to  a  fixed  particle  in  the  fluid. 

Combining  (31e),  (31f),  and  using  &  HR  =  T  £s^,  gives 


9(6T)  _  1  8(6P)  AHR  r 


EXPRESSION  FOR  THE  NET  RATE  OF  REACTION 


(31g) 


For  a  reversible  gas  reaction,  a  rate  equation  of  the  following  type  can  be 


postualted: 


where 


r  -  1 

1  -t?TT 


pyi 

<-A) 

m. 

l 


<  -\  .> 
l 


-T^TT 


Pyi 

(—4 

m. 

i 


<Yi> 


(32) 


<\> 


\  <  0 
\>  0 


T^(T)  and  ^(T)  are  forward  and  backward  rate  constants,  respectively, 
and  as  indicated,  they  are  functions  of  temperature.  The  first  product  term  contains 
the  product  of  the  molal  concentrations  of  the  reactants  raised  to  the  power  of  their 
stoichiometric  coefficients  and  the  second  term  contains  the  products  of  the  reaction. 
At  equilibrium,  F  =  0»  and  so 


and 


1 

Tf(T) 


<  ”Y  i> 


i=l  m. 

l 


_1 _ 

Tb(f) 


—  p  y.  <  y.> 

ii  i  i 


i=l  m. 

l 


(33) 


T^T) 

~f(T) 


(34) 
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KJT)  is  the  equilibrium  constant  of  the  reaction  in  terms  of  molal  concen¬ 


trations. 


By  differentiating  Eq.  (32),  considering  perturbation  terms  to  be  differen¬ 
tials,  and  combining  with  Eqs.  (33)  and  (34),  gives 


\  y 

r _  /  Yl 

.  dlnK  *r . 

r  =  -4 - -  6  T  -  6p 

T  dT  p 


J 

v. 

1=1  1 1  J 


(6F  =  F  since  r  =  0). 

j*/ 

T  ' ,  which  was  defined  by  Eq.  (33),  contains  the  units  of  time  and  must  be 

»i» 

related  to  the  time  of  reaction.  For  a  "very  fast”  reaction,  t  -♦  0,  If  it  is  now 
assumed  that  the  composition  of  the  gas  is  a  continuously  differentiable  function  of 

jJ*  J* 

time,  even  if  t  -*0,  Eq.  (31b)  indicates  that  r  must  remain  finite  as  t  '  -*0.  This 
would  then  imply  that 


j,  r,  [_  Yi  Z 

dlnK  r~.  v 

_^6X  6p  -^6^  =  0 

i=l  1 


for  a  "very  fast"  reaction;  but  this  is  equivalent  to 


f  K  (T)  ]  K  (X) 

-  IT  (Hi>  vil  Tr  Fu  Yi 

' m .  '  'm. 


=  constant. 


The  constant  can  be  evaluated  from  the  situation  at  equilibrium,  Eq.  (34),  and  it 
is  seen  that  the  constant  must  equal  one,  and  so 


z  py. 

"[J  (_i)  Yi  =  K  (T) 
.■  \  'm.  c 

i=l  i 


(38) 
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Therefore,  in  a  ’’very  fast'*  reaction,  the  composition  of  the  gas  mixture  obeys  the 
equilibrium  equation  at  every  instant. 

Eq.  s  (31a),  (31b),  (31c),  (31d),  (31g),  and  (35)  can  now  be  combined  to 
eliminate  fiT,  6p,  and  6y.,  which  gives  three  first  order  linear  partial  differen¬ 
tial  equations  in  three  unknowns. 

These  equations  are 


dink 


(  dT  )  8(6P) 


p  c  at 

M  p  r 


T*ir+\  a^ 
at  A.  1  3xr 

r  i=l  r 


z  .  2 


y.  m.  dlnK 


V  Y-  m. 


ah 


i=i  yi 


dT 


R 


(39a) 


d_p_lP)  8l6P)  + 8(6q)  =i 

p  C  P  T  9t  9x 

K  p  r  r 


m  C 


la 

l  - 


AH 


R 


(39b) 


i_  ®(6P>  +  =  0 


9x 


9t 


(39c) 


EXPLICIT  EXPRESSIONS  FOR  a  and  ar 

e  f 


At  this  point  it  will  be  advantageous  to  determine  explicitly  the  two  sound 
speeds  a^  and  a^. 

Rewriting  Eqs.  (31d)  and  (31e)  in  differential,  rather  than  perturbation 


notation. 
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dP  d0  ,  dT  ,  — 
_  =  —s-  +  —  +  m 


z  j 

V  dyj 

L  ^ 


(40a) 


\r  L.  1 

=  )  S.dy.  +  _£  dT  -  dP 
1  1  T  p  T 


(40b) 


Letting  dS  =  dy.  =  0,  and  then  eliminating  dT  between  Eqs.  (40a)  and  (40b)  gives 


,  3  _  dP  /9p\  _  0 

f  dp  dS  =  0  9p  S,  y.  T  C  "t  P 

J  dy.  =  0  1  P  °  T  “  P 

Ji 


For  gases  obeying  the  ideal  gas  law,  C  -  C  =  ,  where  C  is  the 

Gp  V-  FRoT  V 

constant  volume  specific  heat.  Defining  k  =  and  P  =  - — ,  Eq.  (41)  can  be 

'“’v  m 

put  in  the  form, 


a/2  =  T  —  T 


which  corresponds  in  form  to  the  expression  for  the  speed  of  sound  of  a  nonreacting 
gas.  However,  the  molecular  weight  of  a  nonreacting  gas  remains  constant,  whereas 
for  a  reacting  gas,  it  is  a  variable- 

The  equilibrium  speed  of  sound  can  be  calculated  from  Eqs.  (31b)  and  (36), 
which  are  rewritten  below  in  differential  notation,  a-d  Eqs.  (40a)  and  (40b) 


9(6^)  Y.m. 


(43a) 


^•Y‘  \  Y  'dy- 
i__  do  +  y  -i-i- 


1  >i 


d  InK 

=  ( - -)dT 

dT 


(43b) 


19 


Setting  dS  =  0,  there  is  obtained  by  the  straightforward  combination  of  Eqs,  (40a), 
(40b),  (43a),  and  (43b). 


Z  T  z  Y  ?  m.  AHn  z 

P  Y  i _ L  +  _ RY 


Hd  z  dlnK  T  AH0 

Y-  +  ( _ -> _ - 

P  L  1  dT  7 


M  >  a 


_ _  tT  y  _ 

f  A>i+  TAHR  ,dlnKc,  .  _J_  ( dlnKC' 

'  Z,  w  T5  /IT  -  HT 


I  yi 


The  equilibrium  constant  for  the  reaction  in  terms  of  partial  pressures,  K  , 

P 


is  given  by 


z  pyiR0T  v 

VT)  =  TT.  * 

1  =  1  1 


py.R  T 
o 


The  term -  is  the  expression  for  the  partial  pressure  of  a  perfect  gas  in  a 


gas  mixture. 


K  (T)  and  K  (T)  are  thus  related  by 
c  p  * 


K  =  (R  T) 
c  o  i 


-IY‘ 


dlnK  dlnK 
g  -  _ P 


The  van’t  Hoff  Equation  states  that 


d  InK  AHd 

_ P  _  R 


R  T 
o 
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dlnK 


R  T 
o 


Putting  this  expression  into  Eq.  (44)  and  then  simplifying,  the  equation  for  a  can 
be  put  in  the  form 


■CUT2  z  Y-2m. 

. R.S—  )  -4—1 

m  \r 


“  y. 

1=1  i 


£ Y-a  m. 


£KR2  /r  \2 

+  - —  -  ~C  R  Tf  )  Y  \ 

fn  P  o  I  L  il 

o  ,  ,  ^R2 


T)  — i — -  (C - — )  +  zz - 2  AHL  )  y.+  R  T/)  Y. 

i4l  Yi  Pm  R0T  1  °  \1  1 


In  practice,  the  terms  containing  ^  Y  .  are  negligible  and 

1 


C  RT2r  Y?  m.  A  H2 


pKo  \  Yi  mi 

“  1  7i 


K  v  Y  2  m.  AIL, 

<e„~4  }  >  -4 

P  “  Z-,  “  R  T 


1  'i 


R  T 
o 


If  these  inequalities  are  very  large  so  that  the  terms  containing  AHR2  are 

z  K 

Y  2m 

also  negligible  in  comparison  with  the  )  _i _ i  terms,  Eq,  (50)  becomes 


R  __ 

a  2  =  £  T  =  a  2  . 


v  Y  2m 


The  ratio  of  the  AH^2  term  to  the 


v  t .  in. 

/  — — -  term  in  the  numerator  of  Eq,  (50) 

l  ^i 
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J_  AV 

v  Z  p 

.  r-  Y.  m. 

p  Vs 


The  ratio  in  the  denominator  is 


p  f?i 


y.2m. 

L  1TZ\  _2 _ L 

o  L  Z 


Since,  clearly  C  >  C  -  — -  =  C  ,  the  numerator  ratio  is  less  than  the  denominator 

p  p  v 

m 

ratio  and  so  a  <  a.. 

e  f 


CHARACTERISTICS  OF  THE  GENERAL  SET  OF  DIFFERENTIAL  EQUATIONS 


Having  determined  explicitly  the  expressions  for  a.  and  a  ,  the  character- 

i  e 

istics  of  the  set  of  Eq.  (39)  is  now  calculated,, 

Suppose  5P,  fj  and  §q  are  prescribed  on  a  curve  given  parametrically  by 
xr  =  xr(\)  and  t  =  tr(\).  On  this  curve,  5P  =  6P(M,  F  =  F(M»  6q  =  6q(M.  Then, 
by  the  chain  rule. 


dtr  9UP)  ,  9(SP)  _  d(6P) 

nr  rjt  nr  vz  d\ 

r  r 


(54a) 


dtr  ar  ,  ^r  ar  dr 

nr  w  +  n —  zz  -  zz 

r  r 


(54b) 


dtr  9(&q)  ,  9(§q)  _  d(  §q) 

dX.  5't  dX.  9x  “  dK 

r  r 


(54c) 


Treating  Eqs.  (39)  and  (54)  as  a  system  of  six  linear  equations  in  six  unknowns, 

a(6P)  a(6P)  ar  ar  a(6q)  „  ^  a(sq)  .u  ^ 

— -,  -rr~ — ->  Trr-»  — ,  ,  and  it'---,  the  characteristic  surfaces  of  thes< 


oit 

r 


W’  T’  T5x - 1 

r  9x  r  r 

r 


equations  can  be  found  by  setting  the  determinant  of  the  system,  A,  equal  to  zero. 
Upon  expanding,  the  determinant  is  found  to  be 


*  _  dx  dt  p  C  T-P  dx  _ 

=  •  [(-ar1  <-af>2  -  tL-ibril  • 

p  PCpT 


<lr 


•gsr^i  v  vvw  \  . 


*■■‘"***“55 
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Setting  A  =  0,  it  is  seen  that  if  t  4  0,  characteristics  occur  when 


a) 


dx^  dxr/d\ 

W~  ~  dt  /a\ 
r  r 


=  0,  i.  e. ,  when 


and 


xr  =  x  -  qt 


is  a  constant 


dx  j,  P  C  T 

b) 


i.  e. ,  when 


x  -  (q  ±  a^)t  is  a  constant- 


Smce  the  characteristic  surfaces  propagate  at  the  same  speed  as  the  wave  front, 
condition  (b)  implies  that  the  wave  front  moves  at  speed  af  if  t*  4  0. 


The  compatibility  equation  for  is 


3(6P)  _ 

w~ 


=  -  T 


z 

?• 


'/  \3  AHR"m  CJI 

,&■)( -V-M 


r  + 


_L  (fi,  (f^HdUPJ  +  (^r  sd(6q) 

p*  1  w  +  (-ht  “Hr 


dx 


Compatibility  always  exists  at  =  0.  For  compatibility  at 

r 


(56) 


dx 

r 

ctir 


±  a^  -  ± 


P  C  T 
P 

pC  T-P 
P  P 


/AH  -m  C  T  )  Y. 

»  /  P  V  \„  1  d(6P)  ,  d(6q)  _  „ 

„  r  t  i 1  -  “ar-  ±af-at —  °- 

P  ^p1  /  P  r  r 


(57) 
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9(6y-)  Y-m 

Using  Eq.  (31b),  i.  e.  ,  -  =  — —  T, 

T*  P 


gives 


af2p  d(8yi)  +  1  d(tP)  ±  a  d(8q)  _  Q> 

V  .m.  dt  ^  dt  ^  dt 

2.  i  r  H  r  r 


(58) 


Multiplying  Eq.  (58)  by  p  dt^  and  integrating,  the  compatibility  equation  becomes 


Y.m. 


l  l 


§P  +  a^p  gq  =  constant. 


(59) 


If  t  ‘  =  0,  it  is  seen  that  both  the  determinant  and  the  compatibility  equation 
are  identically  equal  to  zero  and  so  the  wave  front  does  not  necessarily  move  at  the 
speed  a^.  Examination  of  the  system  of  Eqs.  (39)  shows  that  the  highest  derivative 

of  T  is  the  first  derivative  if  t*  f  0.  However,  if  t*  is  equal  to  zero,  the  first 
derivative  term  drops  out  and  the  highest  derivative  of  T  becomes  the  "zeroth" 
derivative.  Thus,  if  T*  is  equal  to  zero,  the  form  of  the  system  of  equations  is 
changed  which  can,  and  in  this  case,  does  produce  a  discontinuous  change  in  the 
characteristics. 

The  characteristics  for  the  case  of  a  "very  fast"  reaction  can  be  found  by 
setting  t  equal  to  zero  in  Eq.  (39a)  eliminating  r  between  Eq.  (39a)  and  (39b)  and 
then  following  a  procedure  similar  to  the  above.  It  is  found  that  in  the  equilibrium 
case,  the  characteristic  surfaces,  or  the  wave  front,  propagate  at  a  speed  ae# 

A 

Since  it  is  physically  impossible  for  t‘  to  ever  actually  be  zero,  these 
results  imply  that  the  wave  front  must  always  propagate  at  a  speed  af.  However, 
the  question  of  the  physical  significance  of  ag  in  near  equilibrium  flow  (t*  very 
small)  is  still  unresolved. 
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DERIVATION  OF  GENERALIZED  ACOUSTIC  EQUATION 

Clarification  of  the  roles  played  by  a^  and  a^  can  be  obtained  by  reducing 
Eq,  (39)  to  a  single  higher  order  equation  and  attempting  to  find  solutions  for  this 
equation.  The  reduction  is  accomplished  in  the  following  way; 

Solve  for  r  in  Eq.  (39b)  to  obtain 


(60) 


Putting  this  Eq.  (60)  for  r  into  Eq.  (39a),  using  Eq.  (41)  and  simplifying,  gives 


If  it  is  assumed  that  6q  and  §P  have  continous  first  derivatives,  then 

=  gljUS}  and  .  Substituting  Eq.  (39c)  into  the  above 

at  dx  0x  9t  ot  ox  ox  9t 

rr  rr  rr  rr 
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equation  gives 


32(6q) 


r  r 


3s(6q) 


r  r 


1  9  UP) 

p  r 


(62a) 


32(6P)  _  32(6P)  _ 


r  r 


r  r 


-  32(6q) 

‘  P  — 

r 


(62b) 


Now,  by  differentiating  Eq.  (61)  with  respect  to  t  ,  combining  with  Eq.  (62a) 

r 

and  simplifying,  using  the  expression  (44)  obtained  for  ae»  there  is  obtained 


1_  92(6P)  _  3gt5P)  +  T**  JLrJL  g^UP)  .  £UP h  =  0 

o  *v  O  O  »  _  L  5  _  p  ^  O  * 


a  2  St  2 
e  r 


at  a-8  at  2  8x  2 
r  f  r  r 


where 


t  =  t 


Z  T 

P 


C  T  y.  m.  TAHd  dlnK  AH 


T-  Z 


7T  T 


~E_  )  _1 — 1+ - 

p  1  yi  p 


p  i  n  \  W 


By  differentiating  Eq.  (61)  with  respect  to  x^,  combining  with  Eq,  (62b)  and 
simplifying,  gives 


i  9zC6q)  3g(6q)  +  T**  8_rJ_  ^(tq)  32(6q) 0_  (65) 

a  3  9t  3  3x  2  3t  af3  3t  3  8x 

e  r  r  r  f  r  r 


2 
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In  x-t  coordinates,  Eqs,  (63)  and  (65)  become 

(q3  X)£_UP)  2q  8a(6P)  ,  1  ^(SP)  + 

a  a  9x 3  a  3  3x8t  a  3  9t 3 

e  e  e 


1(£UP)  +2q  83(6P)  +  1  83(6P)j  =  „ 

ax  a^ 2  3x3t  a^8  3t  3 


(66) 


,  q2  -1}  a3(6q)  +  2q_  33 ( 6q)  +  1  32(6q)  + 

a  2  3x2  a 2  3x3t  a 2  3t2 


'•w*'  , —  O 
t  (q  ^ 


D  83(6q)  f  2q  83(6q)  1  a3(6g)i  =  0 

3x2  a^2  3x3t  a^2  dtz 


(67) 


Eqs.  (63)  and  (65),  or  (66)  and  (67),  can  now  be  regarded  as  generalized  acoustic 
equations  for  the  system. 


PHYSICAL  SIGNIFICANCE  OF  a,  and  a 

f  e 


It  is  seen  that  Eq.  (63)  is  made  up  of  two  wave  equations.  One  wave  equation 

has  characteristics  at  x  ±  a  t  =  constant  or  x  -(q  ±  a  )t  =  constant,  while  the 

r  6  r  6 

other  has  characteristics  at  ±  a^t^  =  constant  or  x  -(q  ±  a^)t  =  constant.  Thus  , 
if  T  '  '  is  exactly  equal  to  zero,  Eq,  (63)  reduces  to  the  equation  of  a  wave  with 

»•*  4>,  »•/  J, 

constant  velocity  a  .  If  t  "is  very  large  (t  -*  oo),  as  would  be  the  case  for  a  very 
slow  reaction  or  for  a  nonreactive  medium,  Eq,  (63)  reduces  to  the  equation  of  a 
wave  with  velocity  a^.  In  any  other  case,  the  equation  is  not  a  wave  equation,  but 
a  third  order  equation.  If  t  is  very  small,  but  not  zero,  the  characteristics  are 
governed  by  the  third  order  terms,  although  the  value  of  the  differential  equation  is 
governed  by  the  second  order  derivatives.  As  has  been  shown,  for  t  ^  0,  no 
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matter  how  small,  the  wave  front  propagates  at  a^.  This  agrees  with  physical 
intuition  in  that  the  wave  front  cannot  have  any  knowledge  of  the  medium  into  which 
it  is  approaching. 

The  significance  of  the  equilibrium  speed  of  sound  will  now  be  determined 
by  attempting  to  find  solutions  of  Eq.  (63).  Since  Eq.  (63)  is  a  linear  differential 
equation  with  constant  coefficients,  one  possible  method  is  to  assume  a  solution  of 
the  form 


oo 


6P 


■I 


CJ 

i 


a.t  +0.(a.)x 
i  r  l  i  r 


i=l 


Here  C.  and  cl  are  undetermined  coefficients  and  (3.  is  a  function  of  a..  By 
putting  this  expression  into  Eq.  (63),  the  relationship  between  (3.,  and  a.  will  be 
obtained.  Then  a.  and  can  be  chosen  as  functions  of  the  integer  i  in  such  a 
way  as  to  satisfy  the  boundary  and  initial  conditions  of  the  particular  problem,  (Note 
that  sinusoidal  solutions  can  be  obtained  if  cl,  |3.,  and  Ch  are  complex). 

Inserting  the  above  expression  into  Eq.  (63)  gives 


a.(3.: 

ri 


!>] 


C.l 

X 


a.t  +p.x 
i  r  r 


=  0. 


(68) 


Requiring  that  each  term  in  the  summation  be  identically  equal  to  zero  implies 


(1  +  cut 


+ 


T 


(69) 
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Note  that  if  **  =  0,  then  -J-  =  ±  a  ,  and  if  t’  *  -*  oo  then  =  ±  a.. 

pT  e’  5  (37  f 

i  ri 

v  v  v'i'3 

If  t  is  small  enough  so  that  t  is  negligible  in  comparison  with  t  , 

*»'  O' 

then  neglecting  second  powers  of  t  , 


❖ 


a.  T  '  a.  a  “ 

ir-*ae[»— i-o-i-g 

1  ^  3.  - 


±  a 


(70) 


a. 

l 


When  t  =  0,  or  t  -*  00 ,  -p—  were  seen  to  be  equal  (both  numerically 
and  dimensionally)  to  the  equilibrium  or  frozen  wave  velocities,  respectively,  which 


were  constants.  If  the  interpretation  of 


K 


as  a  wave  velocity  is  now  extended  to 


the  case  of  t  ,  finite  and  unequal  to  zero,  it  is  seen  that  for  small  t  ‘ ,  the 
"wave  velocities"  are  not  constants  but  weak  functions  of  cu  oscillating  about  the 
value  a  .  Thus,  for  near -equilibrium  flow,  the  expression  of  fiP  can  be  inter- 
preted  as  the  sum  of  a  number  of  disturbances  whose  speeds  are  not,  in  general, 
equal,  but  are  all  close  to  ag.  This  implies  that,  although  the  wave  front  is  mov¬ 
ing  at  a  speed  a^,  the  bulk  of  the  disturbance  behind  the  wave  front  is  moving  at  a 
speed  close  to  a  .  In  order  for  this  to  be  true,  a  must  be  less  than  a.  and  in 
deriving  the  explicit  expressions  for  a^  and  a^,  this  was  shown  to  be  the  case. 

The  significance  of  the  two  sound  speeds,  ag  and  a^,  has  now  been  deter¬ 
mined  for  a  reacting  gas  mixture  near  equilibrium  with  a  fast  reaction  rate.  In 
such  a  system,  a  represents  the  speed  of  the  wave  front,  while  a  represents 
the  approximate  speed  of  the  bulk  of  the  wave  behind  the  front.  Velocity-of -sound 
measurements  have  been  made  in  methane -air  combustion  products  in  which  the 
frozen  speed  a^  was  obtained.  In  some  cases  speeds  greater  than  a^  were 
measured,  and  thes  ^  were  assumed  to  be  close  to  molecular  relaxation. 
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For  the  general  case  of  a  reacting  gas  mixture  near  equilibrium  solutions 
can  theoretically  be  built  up  from  Eq.  (63),  either  by  the  method  shown  or  by  any 
other  technique  for  solving  linear  partial  differential  equations  with  constant  coeffi¬ 
cients.  Practically,  however,  the  usefulness  of  this  procedure  is  limited  by  the 
fact  that  very  little  data  exists  on  rates  of  reaction  and  so,  in  most  cases,  there 

v  V  V 

will  be  no  way  of  determining  and  T£,  and  hence  T  and  t 

A  «A» 

However,  if  t  is  known  to  be  small,  the  bulk  of  the  wave  behind  the  front 
moves  at  an  overall  speed  very  close  to  a^  and  the  acoustic  equation  (63)  is  approxi¬ 
mately  equal  to  the  standard  one -dimensional  wave  equation  with  velocity  a.^.  Thus, 
for  purposes  of  computation,  a  theoretically  instantaneous  reaction  can  be  used  to 
approximate  the  actual  state. 

THE  INSTANTANEOUS  REACTION 

In  an  instantaneous  reaction  the  time  of  reaction  is  zero.  Therefore,  if  the 
gas  flow  is  steady,  the  gas  mixture  will  reach  equilibrium  after  passing  through  an 
infinitesimally  short  distance.  If  a  small  disturbance  is  sent  through  the  pipe,  the 
old  equilibrium  at  any  point  is  destroyed  and  a  new  equilibrium  is  established  at 
exactly  the  instant  that  the  pulse  reaches  that  point. 

Assume  that  the  pulse  starts  at  time  t=0  and  travels  down  the  pipe  with  a 
velocity  a  relative  to  the  gas  stream  or  with  an  absolute  velocity  of  a  +  q.  At  a 
time  t  the  pulse  has  traversed  a  distance  equal  to  (ag  +  q  )t.  For  points 
x  >  (ag  +  q)t,  the  pulse  has  not  yet  arrived  and  the  conditions  are  still  the  old  equi¬ 
librium  conditions.  For  x  <(a  +  q)t,  the  pulse  has  already  passed  and  changed 
conditions  to  the  new  equilibrium.  For  example,  P  is  changed  to  P  +  6P>  where 
6P  is  the  change  in  pressure  between  the  old  equilibrium  and  the  new  Jnd  it  is  of 
perturbation  order.  Similarly,  fcr  x  <(ae+  q)t,  T  =  T  +6T,  q  =  q+6cl> 
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P=  P  +6p,  yi  =  Yi  +  6y^  »  etc. 

According  to  our  hypothesis  of  an  instantaneous  reaction,  for  x  >  (a^+  q)t, 

F=0 since  the  pulse  has  not  yet  arrived,  and  for  x  <(a-e+q)t,  r  =  0  since  the  reaction 
is  already  completed.  At  x  =  (a^  +  q)t,  a  change  of  perturbation  order  is  accom¬ 
plished  in  a  time  interval  approaching  zero.  Hence,  at  x  =  (a^  +  qjt,  rbecomes 
infinite. 

One  is  therefore  led  to  assume  solutions  of  the  form 


6P=6P  U[ae  +  q)t-x]  =?U(agtr-xr) 

(71a) 

6p  =  6p  U  (aetr-xr) 

(71b) 

6y.  =  byAJ  (agtr-xr) 

(71c) 

5q  =  6q  u  (aetr-xr) 

(71d) 

6T  =  fiT  U  (aetr-xr) 

(71e) 

r  =r  Ka  t  -xj 

(71£) 

where  r,  6P,  "sp",  "6y\,  "§q,  6T  are  constant  amplitude  terms.  "U"  is  a  symbol 
for  the  unit  (Heaviside)  step  function  defined  by  U(\)  =  1  for  X  >  0  and^UW  =  0 
for  X  <  0.  The  impulse  function  is  defined  by  I(X)  =  0,  if  X  ^  0  and  I(X)dX  =  1. 
Although  neither  U(X)  nor  I(X)  is  an  everywhere  continuous  function 
and  I(X)  is  not  properly  a  function  at  all  ,  they  can  nevertheless  be  treated  as 
continuous  differentiable  functions  which  are  related  by  I(\)  —  1  • 

Putting  Eq.  (71)  into  Eqs.  (31a),  (31b),  (31c),  (31d),  (31g),  and  (36),  which 
are  the  governing  equations  for  an  instantaneous  reaction,  gives 

ae‘6p"  -  o  6p  =  0 


(72a) 
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7^7  6Yi  -  T  =  °  (72b) 

i  i 


p  ag  6q  -  6P  =  0 


(72c) 


iE.  +  +35  \ 

P  T  L 

1 


(72d) 


ae6T 


_  AH 

6P 

P  c 


R 

P 


0 


d’ln'K  _ 

( - -)6T  - 

dT 


P 


(72e) 


(72f) 


Eqs.  (72)  are  a  set  of  linear  homogeneous  equations  in  six  unknowns.  The 
requirement  for  the  existence  of  non-trivial  solutions  is  that  the  determinant  of  the 

system  be  equal  to  zero  and  evaluation  of  the  determinant  shows  that  this  is  true  if 
ag  is  given  by  Eq.  (44).  This  being  the  case,  Eqs.  (72)  enable  one  to  obtain  all  the 
perturbation  amplitudes,  if  any  one  of  them  is  known. 
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